Cosmological simulations using a static scalar-tensor theory 
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We present ACDM iV-body cosmological simulations in the framework of a static general scalar- 
tensor theory of gravity. Due to the influence of the non-minimally coupled scalar field, the gravita- 
tional potential is modified by a Yukawa type term, yielding a new structure formation dynamics. 
We present some preliminary results and, in particular, we compute the density and velocity profiles 
of the most massive group. 
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I. INTRODUCTION 

The problem of explaining the structure formation of 
the Universe is one of the most fascinating at the be- 
ginning of this new millennium. From the recent exper- 
imental developments and observations, cosmology has 
acquired the status of a high precision area of research 
In fact, recent and independent observational data 
measured in the CMBR on various angular scales @, 
type la supernovae [4] , as well as in galaxy surveys 
suggest that ft = Q\ + Fl m « 1, or that £7 a « 0.73 and 
Q m « 0.27, implying the existence of dark energy and 
dark matter, respectively. However, the nature of dark 
components is still unknown. These observations favour 
the now standard ACDM cosmological model. Despite 
some problems on galactic scales revealed by numerical 
simulations, which could be explained with some alterna- 
tive models (e.g. Warm Dark Matter), large scale struc- 
ture simulations, predicted by the ACDM model, agree 
well with observations. Naturally, particular inflationary 
scenarios motivated by different particle physics theories 
have their own dark matter candidates, such as axion, 
neutralino, Higgs particle, among others, as well as a 
quintessence field @, Q . In order to include such parti- 
cles in cosmological models one usually adds a scalar field 
(SF) equation to the general relativity field equations. 
One possibility is to couple this field non-minimally to 
gravity to have a scalar-tensor theory of gravity (STT). 

In this work we present some results of the role played 
by a massive non-minimally coupled SF on the ACDM 
universe structure formation process. Our theoretical 
model is built from a general STT static SF which mod- 
ifies the cosmological dynamics and the Newtonian grav- 
ity potential on particles; the dynamical full treatment of 
SF perturbations is now under development. We evolve a 
cosmological cube using standard ACDM equations with 
periodic boundary conditions, where the particles inter- 
act through the Newtonian force plus an additional term. 
The latter comes from a Yukawa type potential derived 
from the Newtonian limit of a STT. It should be noted 
that in the present work the SF does not replace the dark 



matter or dark energy, but rather coexists with them. To 
perform the simulations we have modified a standard se- 
rial treecode developed by one of us (MARM) @ and 
Gadget 1 [10] in order to take into account the contribu- 
tion of the Yukawa potential. 



II. EVOLUTION EQUATIONS USING A STATIC 

STT 

In a previous paper we found the solutions for the 
potential-density pair problem in the Newtonian limit of 
a scalar-tensor theory of gravity [ll|. Here, we applied 
those results and found that the potential of a single par- 
ticle of mass m is given by [l2[ 
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For local scales, r <C A, deviations from the Newtonian 
theory are exponentially suppressed, and for r ^> A the 
Newtonian constant diminishes (augments) to Gat/(1 + 
a) for positive (negative) a. This means that equation 
(PQ) fulfills all local tests of the Newtonian dynamics, and 
it is only constrained by experiments or tests on scales 
larger than -or of the order of- A, which in our case is of 
the order of galactic scales. 

To simulate cosmological systems, the expansion of the 
Universe has to be taken into account. Here, we employ a 
cosmological model with a static scalar field which is con- 
sistent with the Newtonian limit given by Eq. (pQ). Thus, 
the scale factor, a(t), is given by the following Friedman 
model, 



a 6 H 2 



Hi 



a 



mO 



^ao a 3 



1 



a 



m0 



AO 



1 



a 



(2) 

where H = d/a, flmo and Q\o are the matter and energy 
density evaluated at present, respectively. We notice that 
the source of the cosmic evolution is deviated by the term 
1 + a when compared to the standard Friedman-Lemaitre 
model. Therefore, it is convenient to define a new den- 
sity parameter by Qf = ^/(l + a). This new density 
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parameter is such that + = 1, which implies a 
fiat Universe, and this shall be assumed in our following 
computations, where we consider (0^,0^) = (0.3,0.7). 
For positive values of a, a flat cosmological model de- 
mands to have a factor (1 + a) more energetic content 
(O m and ft \) than in standard cosmology. On the other 
hand, for negative values of a one needs a factor (1 + a) 
less Qm and £7 a to have a flat Universe. To be consis- 
tent with the CMB spectrum and structure formation 
numerical experiments, cosmological constraints must be 
applied on a in order for it to be within the range (—1,1) 
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In the Newtonian limit of STT of gravity, the Newto- 
nian motion equation for a particle i is written as 

^ +2H ^ = £ it^f ^(i*-x,i,«, a) 
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(3) 

where x is the comovil coordinate, and the sum includes 
all periodic images of particle j, and Fsf(t, a, A) is 

F SF (r, a, A) = 1 + a (l + e~ r ' x (4) 

which, for small distances compared to A, is Fsf(t < 
A, a, A) ~ 1 + ol (l + jj and, for long distances, is 
Fsf(^ > A, a, A) « 1, as in Newtonian physics. 

We now analyze the general effect that the constant 
a has on the dynamics. The role of a in our approach 
is as follows. On one hand, to construct a flat model 
we have set the condition fi^ + = 1, which implies 
having (1 + a) times the energetic content of the stan- 
dard ACDM model. This essentially means that we have 
an increment by a factor of (1 + a) times the amount 
of matter, for positive values of a, or a reduction of the 
same factor for negative values of a. Increasing or re- 
ducing this amount of matter affects the matter term on 
the r.h.s. of the equation of motion (j3]), but the amount 
affected cancels out with the term (1 + a) in the denomi- 
nator of ([3|) stemming from the new Newtonian potential 
(PQ). On the other hand, the factor Fsf augments (di- 
minishes) for positive (negative) values of a for small dis- 
tances compared to A, resulting in more (less) structure 
formation for positive (negative) values of a compared to 
the ACDM model. For r ^> A the dynamics is essentially 
Newtonian. 



III. RESULTS 

We now present results for the ACDM model of the 
Universe model previously described. Because the visi- 
ble component is the smaller one and given our interest 
to test the consequences of including a SF contribution 
to the evolution equations, our model excludes gas par- 
ticles, but all its mass has been added to the dark mat- 
ter. Therefore, our model is as follows. We start our 
simulation with an initial distribution of N = 2 x 32 3 
particles in a box with sides of 50/z -1 Mpc at z = 10. 



This case is similar to the one that comes with Gadget 1 
H3|. At present epoch, fi£ = 0, fi^ = 0.3, Q% = 0.7, 
Ho = lOOh km/s/Mpc, h = 0.7. We restrict the values 
of a to the interval (-1, 1) 0, 0J, [H, and use A = 5 
Mpc, since this scale turns out to be an intermediate scale 
between the size of the clump groups and the separation 
of the formed groups. 

In Fig. 1 we show y-z snapshots at redshift z = of 
our ACDM model. Fig. 1 (a) presents the standard case 
without SF, i.e., the interaction between bodies through 
the standard Newtonian potential. In (b) we show the 
case with a = 1, A = 5 Mpc. In (c) a = — 1/2, A = 5 
Mpc. In (d) a = —1/4, A = 5 Mpc. One notes clearly 
how the SF modifies the matter structure of the system. 
The most dramatic cases are (b) and (c) where we have 
used a = 1 and a = —1/2, respectively. Given the ar- 
gument at the end of last section, in the case of (b), 
for r <C A, the effective gravitational pull has been aug- 
mented by a factor of 2, in contrast to case (c) where 
it has diminished by a factor of 1/2; in model (d) the 
pull diminishes only by a factor of 3/4. That is why one 
observes for r < A more structure formation in (b), less 
in (d), and lesser in model (c). The effect is then, for a 
growing positive a, to speed up the growth of perturba- 
tions, then of halos and then of clusters, whereas negative 
a values (a — > —1) tend to slow down the growth. 

Next, we found the groups in the system using a friend- 
of-friend algorithm and select one of the most massive 
ones. The chosen group is located approximately at y = 
19 Mpc, z = 12 Mpc, and it is labeled with the letter "G" 
in Fig. 1(a). The group was analyzed by obtaining their 
density profiles and circular velocities. In Fig. 2(a) we 
show the density profiles for this group. The more cuspy 
case is for a = 1 and the less cuspy is for a = —1/2. 

In Fig. 2(b) we show, for the same group, circular 
velocity curves, computed using v 2 c = GjsrM(r)/r. The 
case with a = 1 corresponds to higher values of v c , since 
this depends on how much accumulated mass there is at 
a distance r and this is enhanced by the factor Fsf for 
positive values of a. 

The groups are at most 2 Mpc in size at z = 0, which 
means that the inner structure is such that r < A, being 
affected by the factor Fsf as explained at the end of 
section 2. While for r > A the overall structure formation 
process is governed by Newtonian physics, which is the 
reason why the overall structure of the models in Fig. 1 
is similar. 



IV. CONCLUSIONS 

We have used a general, static STT that is compatible 
with local observations by the appropriate definition of 
the background field constant, i.e. < (j) >= G^ 1 ^ + a). 
A direct consequence of our approach is that the amount 
of matter (energy) has to be increased for positive values 
of a and diminished for negative values of a with respect 
to the standard ACDM model in order to have a flat cos- 
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(xlO 4 ) (a) ACDM - z=0 - without SF 
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(xio 4 ) (b) ACDM - z=0 - SF: o=l, X=5 Mpc 
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(xio 4 ) (c) ACDM - z=0 - SF: a=-l/2, X=5 Mpc 
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(xio 4 ) (d) ACDM - z=0 - SF: a=-l/4, X=5 Mpc 
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FIG. 2: (a) Density profiles for one of the most massive groups 
at z = of a ACDM universe. The group is located approxi- 
mately at y = 19 Mpc, z = 12 Mpc, labeled with "G" in Fig 
1(a). Vertical scale is in units of p = lO^Mofc-V^kpc) 3 . 
(b) The corresponding circular velocity. 



mological model. Quantitatively, our model demands to 
have fi/(l + a) = 1 and this changes the amount of dark 
matter and energy of the model for a flat cosmological 
model, as assumed. 

The general gravitational effect is that the interaction 
including the SF changes by a factor FsF(r, a, A) « 1 + 
a (l + jj for r < A in comparison with the Newtonian 
case. Thus, for a > the growth of structures speeds up 
in comparison with the Newtonian case. For the a < 
case the effect is to diminish the formation of structures. 
For r > A the dynamics is essentially Newtonian. 

Using the resulting modified dynamical equations, we 
have studied the structure formation process of a ACDM 
universe. We varied the amplitude and sign of the 
strength of the SF (a) in the interval (-1,1) and per- 
formed several 3D-simulations with the same initial con- 
ditions. From our simulations with different values of a, 
we have found that the inclusion of SF changes local dy- 



FIG. 1: y—z snapshots at z = of a ACDM universe. See text 
for details. 
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namical properties of the most massive group considered, 
and accordingly the density profile and circular velocity, 
however, the overall structure is somewhat similar. Here, 
we notice that we have also studied other massive groups, 
in particular one smaller group located approximately at 
y = 16.5 Mpc, z = 22 Mpc in Fig. 1(a). The trends 
are quite similar to the most massive group, so that our 
conclusions prevail. 

In this work we only varied the amplitude of the SF 
(a) leaving the scale length (A) of the SF unchanged. 
After some preliminary runs, the increase of A enhances 
the structure formation process for a positive, and the 



decrease of A makes the structure grow at a slower rate. 
In a future work we will study an ampler space parameter 
and with a much better resolution. Also, the cosmologi- 
cal initial conditions will be constructed using the matter 
density field corresponding to the modified gravity. 
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